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@) Continuous display of peak and mean blood flow velocities. 



(57) A technique for continuously determining 
and displaying the peak and mean velocities of 
spectral Doppier information is disclosed. Re- 
ceived spectral Doppier data is examined over a 
predetermined period to determine the signal to 
noise ratio of the Doppier signals. Individual 
instantaneous measurements, or lines, of Dop- 
pier spectral information are analyzed in the 
power versus frequency domain to compare the 
line information against a noise threshold deter- 
mined by the signal to noise ratio. The 
maximum frequency signal component is there- 
by identified as the peak velocity of the spectral 
line. Mean velocity is calculated as an intensity 
weighted mean about the maximum power fre- 
quency of the line data. The identified peak and 
mean velocities are displayed in contrasting 
colors or shades on a continuously updated and 
displayed flow velocity versus time spectral dis- 
play. 
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This invention relates to ultrasonic diagnostic 
systems which measure the velocity of blood flow us- 
ing spectral Doppler techniques and, In particular, to 
the continuous display of such information including 
peak and mean blood flow velocities. 5 

Ultrasonic diagnostic imaging systems are in 
widespread use for performing ultrasonic imaging 
and measurements through the use of scanheads 
which are placed against the skin. Such systems are 
in common use by radiologists, cardiologists, and ob- to 
stetri cians for examinations of the heart a developing 
fetus, or internal abdominal organs. These systems 
operate by actuating an acoustic transducer within 
the scanhead to transmit ultrasonic energy through 
the skin and into the body and receiving ultrasonic 15 
echoes returned from the internal structure of the 
body. The transducer converts the acoustic energy of 
the returning ultrasonic echoes into electrical signals 
which are coupled by a cable to the diagnostic system 
for processing and display. 20 

In the measurement of blood flow characteristics, 
waves of returning ultrasonic signals are compared to 
a phase reference to determine the phase shift of the 
returning waves. As the transmitted ultrasonic wave 
impinges upon flowing material such as blood cells, 2s 
the movement of the flowing material will impart a 
Doppler shift to the returning echo signal. This phase 
shift which is commonly measured in kilohertz, trans- 
lates into the rate of movement or velocity of the blood 
flow. This Doppler velocity information is convention- 30 
ally displayed as a continuous spectrum of lines of 
varying amplitudes in a moving or scrolled display on 
a video monitor. Each line represents an instantane- 
ous measurement of blood flow. For arterial flow, the 
spectrum of lines will rise during the systolic phase of 35 
the heart cycle and decline to lower levels during the 
diastolic phase of t he heart cycle. As the flow of blood 
in a vessel or the heart is continuously monitored and 
interrogated by Doppler ultrasound, the spectrum of 
systolic and diastolic velocities is continuously dis- 40 
played and passed before the clinician. 

In order to analyze the various disease states of 
the cardiovascular system it is desirable to calculate 
a number of parameters from the spectral velocity in- 
formation. These parameters include peak systolic 45 
velocity, minimum diastolic velocity, the systolic/dia- 
stoiic ratio, the pulsatility index, the velocity time in- 
tegral, and the time averaged mean velocity. Many of 
these parameters are predicated upon identifying 
peak and mean velocities of blood flow. The conven- so 
tional technique for determining the peak and mean 
velocities and related parameters is to save a record- 
ing of a period of spectral information over several 
consecutive heart cycles. With the spectral display 
held motionless on the screen the clinician manually 55 
traces the peaks of the spectral display with a cursor. 
Calculation software may then use the spectral trac- 
ing to calculate mean velocities and a variety of other 



parameters. 

This technique has several obvious disadvantag- 
es including the need to do a tedious, time consuming 
tracing and the inaccuracies inherent in a manually 
executed procedure. Moreover, this technique does 
not lend itself to automated processing, as the calcu- 
lation software is dependent upon this manual exer- 
cise of display tracing. Furthermore, the ultrasonic 
examination of a patient must be interrupted by stop- 
ping the acquisition of Doppler information in order to 
perform the manual tracing of previously acquired 
data. It would be preferable to be able to automatically 
acquire the peak and mean velocity data without the 
need for manual intervention or interruption of a pa- 
tient examination. 

A number of problems confront any attempt to au- 
tomatically analyze Doppler spectral data, however. 
Any such technique must account for and respond ap- 
propriately to faulty spectral information or loss of 
Doppler signals which may occur when a patient 
moves or the clinician moves the scanhead, for in- 
stance. Also, the technique must be able to reliably 
distinguish between Doppler signal information and 
noise so that the peaks of weak signals are accurately 
determined in the presence of background noise. 

In accordance with the principles of the present 
invention a technique for continuously determining 
and displaying the peak and mean velocities of spec- 
tral Doppler information is disclosed. Received spec- 
tral Doppler data is examined over a predetermined 
period to determine a noise threshold level. Individual 
instantaneous measurements, or lines, of Doppler 
spectral information are analyzed to compare the line 
information against the noise threshold. The correla- 
tion of the threshold and the spectral line information 
identifies the peak velocity of the spectral line. Mean 
velocity is calculated as an intensity weighted mean 
of the spectral line information. The identified peak 
and mean velocities are displayed in contrasting col- 
ors or shades on a continuously updated and dis- 
played frequency versus time spectral display. 

In the drawings: 

FIGURE 1 illustrates a scrolling Doppler spectral 
display as produced by ultrasound systems of the 
prior art; 

FIGURE 2 illustrates a power versus frequency 
plot of the data points of a typical Doppler spec- 
tral measurement 

FIGURE 3 illustrates the identification of the 
peak velocity of a spectral line in accordance with 
the principles of the present invention; 
FIGURE 4 illustrates a continuous display of 
spectral Doppler peak and mean velocities in ac- 
cordance with the principles of the present inven- 
tion; and 

FIGURE 5 illustrates a preferred monochrome 
display of continuous peak and mean velocity in- 
formation. 
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Referring first to FIGURE 1, a scrolling Doppler 
spectral display which is produced by ultrasound sys- 
tems of the prior art is shown. The illustrated display 
is developed by repetitively transmitting ultrasonic 5 
Doppler waves to a selected site, such as the heart 
or a blood vessel, inside the body. Echo signals are 
returned by moving blood cells In the heart or vessel 
and received by a piezoelectric transducer which con- 
verts the ultrasonic echoes into electrical signals. The 10 
signals are amplified and phase detected to deter- 
mine their frequency shift characteristics. Samples of 
the detected signals are processed in a Doppler proc- 
essor to determine the power (intensity) versus fre- 
quency characteristics of the signals. The spectral 15 
frequency characteristics are translated to velocity 
equivalents, and the Doppler information of discrete 
sampling periods is displayed as a sequence of con- 
tinuous scrolling spectral lines in a realtime time ver- 
sus velocity display as shown in FIGURE 1 . In the dis- 20 
play of FIGURE 1 , newly generated spectral lines are 
produced at the right side of the display. The se- 
quence of lines moves or scrolls from right to left, with 
previously generated spectral data on the left and 
progressively more current data to the right Each line 25 
conveys the range of flow velocities detected in the 
blood flow at a particular time of Doppler interroga- 
tion. The highest velocities shown bylines 10, 20, and 
30 would typically occur during the systolic phase of 
the heart cycle. The intervals 1 2, 22, and 32 between 30 
the systolic phases represent flow velocity during the 
intervening diastolic phases of heart action. 

In a typical diagnostic procedure of the prior art 
the clinician manipulates the ultrasonic scan head 
and steers the ultrasonic beam toward the vessel or 35 
organ where flow velocity information is desired. The 
spectral display is monitored as its scrolls by until the 
clinician is satisfied that it has become stable. The 
spectral display is then frozen on the screen and 
saved for analysis. The analysis may proceed by stop- 40 
ping the scanning of the patient and manually tracing 
the spectral peaks with a cursor controlled by a joy- 
stick or trackball on the ultrasound system. Calcula- 
tion software in the system can then operate on the 
tracing to determine clinical flow parameters such as 45 
peak systolic velocity, minimum diastolic velocity, the 
systolic/diastolic ratio, the pulsatility index and the 
velocity time integral. The time averaged mean veloc- 
ity can then be estimated by operating on the peak 
velocity tracing data in concert with assumptions so 
made as to certain flow characteristics. Alternatively 
the saved spectral information can be applied to an 
offline processor which is capable of operating on the 
spectral information to determine desired clinical 
parameters. In either case, patient scanning is inter- 55 
rupted and time consuming operations must be per- 
formed to derive the desired diagnostic information. 

The display of FIGURE 1 may be produced by us- 
ing the spectral Doppler acquisition and display capa- 



bilities of a number of commercially available ultra- 
sound systems such as the Ultramark(R) 9 diagnostic 
ultrasound system manufactured by Advanced Tech- 
nology Laboratories, Inc. of Bothell, Washington. 

FIGURE 2 illustrates a typical intensity versus 
frequency plot of Doppler data produced by the ultra- 
sound system's Doppler processor from signals ac- 
quired during a particular time of Doppler interroga- 
tion of the flow of the vessel or organ being examined, 
in a digital signal processing system such as the Ul- 
tramark 9 the Doppler data is a series of discrete dig- 
ital values as denoted by the circles and Xs of FIG- 
URE 2. For ease of illustration a curve 40 in FIGURE 
2 has been drawn to connect the digital values. The 
curve 40 and its digital values will range about a vert- 
ical axis 42, extending over a range of frequency val- 
ues in the horizontal direction. The vertical axis 42 
marks a frequency value of zero on the horizontal fre- 
quency axis. The frequency limits vary with the sys- 
tem setting for the range of velocities to be detected. 
Exemplary limits might be +1000 Hz and -1000 Hz, 
which correspond to flow velocities of +30 cm/sec 
and -30 cm/sec, respectively. In the vertical intensity 
direction the point P(f)max denotes the maximum pow- 
er or intensity of the received signals at its corre- 
sponding frequency on the horizontal frequency axis. 
The frequency of peak power is not the peak signal 
frequency, however, which it is the objective of the 
present invention to identify. 

To positively identify the peak signal frequency, 
valid Doppler signals must be distinguished from 
noise. Without a viable noise immunity technique, a 
peak signal tracing technique can erroneously identi- 
fy noise peaks as signal peaks. Thus, in accordance 
with the principles of the present invention, a signal 
to noise threshold for the received signals is first de- 
termined. Signal to noise determination begins by se- 
lecting a number of spectra] lines during a heart cycle 
for analysis. While FIGURE 1 shows approximately 
ten spectral lines during one heart cycle, such as the 
ten lines between systolic peaks 10 and 20, in prac- 
tice the number of such lines is generally far greater, 
ranging up to fifty or more. The number of lines is a 
function of several factors such as the PRF (pulse 
repetition frequency) of Doppler interrogation, signal 
averaging, and display scroll rate. From the total num- 
ber of lines occurring during the heart cycle only ev- 
ery third, or every sixth, or every seventh line may be 
selected for analysis. Alternatively, lines may be se- 
lected over some predetermined period, such as ev- 
ery fourth line occurring over a period of one second. 
Depending upon the selection methodology chosen, 
a number of lines such as sixteen or thirty- two may 
be selected during the chosen interval. 

For each of these spectral lines, each having a 
unique data set represented by the one of FIGURE 2, 
the peak intensity P(f) max is found. Then an assumed 
noise threshold is applied to each data set as indicat- 
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ed by the threshold NOISE,™ in FIGURE 2. This as- 
sumed threshold couid be 3 dB below the peak inten- 
sity P(f)max» for instance. Then all data points above 
and below the NOISE max threshold are respectively 5 
averaged for each line. In the example of FIGURE 2, 
the values of the X data points are averaged to deter- 
mine an average signal value for the spectral line and 
the values of the O data points are averaged to de- 
termine an average noise value for the spectral line. 10 
In the determination of the average noise value, val- 
ues of zero may be ignored such as those which 
would occur on either side of the origin in FIGURE 2. 
Next, the average signal values of all of the lines are 
averaged to determine an average signal value for 15 
the selection period and the average noise values of 
all of the lines are averaged to determine an average 
noise value for the selection period. From these se- 
lection period average values a signal to noise ratio 
SNR is calculated as 20 

SNR = Avera 9 eSi 9 na U 
k * Average Noise^ 
where the constant k Is a function of the scroll rate, 
or display rate, of the spectral lines. The calculated 
value of SNR is in units of dB [down from the signal 25 
intensity peak]. 

Thus, at approximately the end of a first heart cy- 
cle of spectral lines, a signal to noise value has been 
determined. This SNR value is used as explained be- 
low to identify the peak velocity of each spectral line 30 
during the succeeding heart cycle. Concurrently a 
new SNR value is calculated in the foregoing manner 
during this second heart cycle for use in identifying 
peak velocities during the third heart cycle. 

During the second heart cycle each spectral line 35 
is analyzed and its peak velocity value identified as 
illustrated with reference to FIGURE 3. This FIGURE 
shows a curve 50 which connects the discrete digital 
values of a spectral line, some of which are indicated 
by the circles on the left side of the curve. The first 40 
step in peak velocity identification is to find the peak 
intensity P(f) max of each successive spectral line and 
analyze the P(f)max value to determine whether the 
spectral line is a valid one. This analysis provides the 
peak signal determination with a degree of noise im- 45 
munity. 

The Pfflmax value is first compared to the value of 
the average noise of the previous heart cycle, Aver- 
age Noise,*. If the value of P(f)max is less than the val- 
ue of the average noise, no peak determination is so 
made for this spectral line; the line is ignored and the 
analysis proceeds to the next spectral line. If five suc- 
cessive spectral lines are analyzed in this manner 
and fail the comparison, the value of Average Noise,* 
is decremented by a predetermined amount and the 55 
analysis continues using the decremented value of 
Average Noise,*. If no valid comparison is found and 
the value of Average Noise,* is decremented below a 



predefined lower limit, the system concludes that 
only noise is being received; no peak velocity deter- 
minations are made until the level is exceeded. Such 
a condition could arise if the scanhead was removed 
from the patient for instance, or if the Doppler beam 
is not being directed to a region of the body where 
flow is occurring. 

If the Pt^max value meets the average noise com- 
parative test, the PtfJnw value is then analyzed to see 
if it falls within a range of possible values. If Pff),^ is 
below the range the spectral line is probably noise in- 
duced and is discarded. If P(f) max is within the range 
the signal is probably bordering on noise and cannot 
reliably be identified as valid spectral data; If this oc- 
curs the previously identified peak value of the previ- 
ous spectral line is held over and analysis proceeds 
to the next spectral line. If the PffW value is above 
the range the signal is strong enough to be regarded 
as valid data for the performance of a new peak ve- 
locity determination. The preferred end limits for the 
range are functionally of the form Kj.* Average Noise,* 
and K H *Average Signal,*, representing the lower and 
higher range limits, respectively. 

Once the P(f) max value of a particular spectral line 
satisfies the foregoing noise analysts and the line has 
therefore been found to exhibit the desired immunity 
to noise and spurious effects, the value of P(f)max is 
examined to see which side of the vertical baseline 52 
the value is located. This determines the direction of 
the flow in relation to the position of the scanhead, 
and also the end of the spectrum from which the 
search for the peak velocity is to proceed. In the ex- 
emplary spectral line of FIGURE 3, the Pif),^ value 
is located to the left of the baseline 52, indicating flow 
toward (or away from, depending upon the chosen 
convention) the scanhead. The search for the peak 
velocity value therefore begins from the leftmost dig- 
ital value, which is circle 54. 

The search for the peak velocity value proceeds 
from value 54 and continues to the right in the draw- 
ing, through value 56 and then to the value indicated 
by the solid circle 58. Between values 56 and 58 the 
SNR threshold is traversed. The system will choose 
the value closest to the SNR threshold as the peak 
velocity value, which in the example shown is digital 
value 58. This value has a frequency of f p as shown 
in the drawing, which is substantially the highest fre- 
quency of the values which are above the SNR 
threshold. The velocity corresponding to fp is thus 
identified as the peak velocity for this spectral line, 
and that velocity is graphically marked in the spectral 
line display. FIGURE 4 illustrates a spectral line dis- 
play in which the peak velocity of each spectral line 
has been identified in this manner and the peaks con- 
nected by the solid display line 60. As FIGURE 4 
shows, the spectral line peak velocities can be iden- 
tified and displayed as the spectral lines occur and 
are displayed, thereby providing a realtime continu- 
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ous display of peak spectral velocities. 

For each displayed spectral line which satisfies 
the noise immunity tests a mean velocity value is cal- 
culated and displayed. A variety of techniques are 5 
known for calculating mean velocity, a preferred one 
being an intensity weighted mean of the spectral P(f) 
digital values. This calculation is of the form 

Mean Velocity = ZK ™ *J P(f) 

1 X P(f) 10 

where K„, is a function of the location of P(f) max on the 
frequency axis. A correction is then made for the ef- 
fects of noise, which may be any of a variety of accept- 
ed correction techniques such as those described in 
the paper by R. W. Gill entitled "Performance of The 15 
Mean Frequency Doppler Modulator", published in Ul- 
trasound In Med. & Biol. , Vol 5, pp. 237-47 (1 979). The 
mean velocity thereby determined is marked on the 
spectral line display, also concurrently with the initial 
appearance of the spectral line at the right-hand side 20 
of the spectral line display. FIGURE 4 shows a dash- 
ed line display which connects the calculated mean 
velocity values of the displayed spectral lines. 

Several preferred techniques may be employed 
for displaying the spectral lines with the peak and 25 
mean velocity determinations. For a color display the 
peak velocity line 60 and the mean velocity line 62 are 
displayed in colors or shades which contrast with the 
spectral lines and with each other. One of the peak 
or mean velocity lines may be emphasized over the 30 
other, the preferred technique being to more brightly 
display the peak velocity line 60 as compared to the 
mean velocity line 62, as by using a greater pixel line 
width. FIGURE 5 illustrates a preferred technique for 
a monochrome display. (NOTE: In FIGURE 5 the 35 
shadings are reversed from what is described.) In this 
FIGURE the spectral lines 70 are displayed in shades 
of grey against a black background 72. The peak ve- 
locity line 80 is brightly displayed in white and the 
mean velocity values are marked by blanking their 40 
positions on the respective spectral lines, thereby ef- 
fectively leaving a black line running through the the 
spectral lines 70 as indicated at 82. This technique 
takes advantage of rapid, high density production and 
display of spectral lines, in which the spectral lines 70 45 
are displayed virtually adjacent to each other, thereby 
resembling a continuous band of grey shading below 
the peak velocity line 80. The black mean velocity line 
82 is thus distinctly displayed in contrast to the sur- 
rounding grey shading of the spectral lines. so 

A number of additional techniques may be em- 
ployed to insure greater integrity for the displayed 
Doppler velocity information and noise immunity. For 
instance, if spectral lines are acquired at a rate in ex- 
cess of the display scroll rate, a number of acquired 55 
lines may be averaged and the average of several sig- 
nal acquisitions displayed as one spectral line. The 
number of excessive changes in P(f) rnax values within 



a given time interval can be monitored. The occur- 
rence of a number of excessive peak changes within 
a given interval indicates that random noise signals 
are being processed, prompting the system to stop 
the display of peak and mean velocities until valid 
spectral line data is identified. The identified peak 
and/or mean velocity values of a number of succes- 
sive spectral lines can be examined and the average 
or median value displayed to further diminish noise 
and other random effects. These further enhance- 
ments will increase the value of the displayed infor- 
mation to the diagnosing clinician. 

Since the technique of the present invention 
presents continuous realtime peak and mean velocity 
information while the clinician is scanning the patient, 
a variety of clinical parameters which require peak or 
mean velocity information for their determination can 
be concurrently calculated and numerically displayed 
for the clinician. These include peak systolic velocity, 
end diastolic velocity, minimum diastolic velocity, 
time averaged peak velocity, systolic/diastolic ratio, 
pulsatility index, velocity time integral, resistance in- 
dex, pressure half time, acceleration time, decelera- 
tion time, heart cycle length, and time averaged mean 
velocity. A listing of such conventional parameters is 
given in the text Doppler Ultrasound by D.H. Evans ef 
a/., pub. by John Wiley & Sons (1989) at pages 166- 
84 and 266-69. With this further capability a diagnos- 
tic examination of cardiovascular performance can 
be performed reliably and efficiently to the benefit of 
both patient and clinician. 



Claims 

1. Ultrasonic diagnostic apparatus for providing 
Doppler information concerning flow velocities 
comprising: 

means for acquiring Doppler signals from 
a region of a body; 

means for processing said Doppler signals 
to produce spectral line data for continuous dis- 
play of a succession of spectral lines during the 
continuation of Doppler signal acquisition; 

means for analyzing the data of a plurality 
of said spectral lines to determine a threshold lev- 
el; 

means for utilizing said threshold level and 
said spectral line data of a first successive spec- 
tral line to identify a peak velocity value in said 
first successive spectral line prior to the display 
of said first successive spectral line; and 

means for displaying said first successive 
spectral line and its identified peak velocity value 
prior to the display of a second successive spec- 
tral line. 

2. The ultrasonic diagnostic apparatus of Daim 1, 
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further comprising: 

means for utilizing said spectral line data 
of said first successive spectral line to identify a 
mean velocity value in said first successive spec- 5 
tral line prior to the display of said first successive 
spectral line; and 

means for displaying said first successive 
spectral line and its identified mean velocity value 
prior to the display of a second successive spec- 1 o 
tral line. 



of grey against a black background, said identi- 
fied peak velocity values brighter than said 
shades of grey, and said mean velocity values in 
said background color. 



3. The ultrasonic diagnostic apparatus of Claim 2, 
wherein said means for utilizing said spectral line 
data of said first successive spectral line to Iden- is 
trfy a mean velocity value comprises means for 
calculating an intensity weighted means utilizing 

said spectral line data of said first successive 
spectral line. 

20 

4. The ultrasonic diagnostic apparatus of Claim 1 , 
wherein said means for analyzing the data of a 
plurality of said spectral lines comprises means 
for determining a noise threshold level. 

25 

5. The ultrasonic diagnostic apparatus of Claim 4, 
wherein said means for analyzing comprises 
means for analyzing the data of a plurality of said 
spectral lines occurring during a first time interval 

to determine a threshold level; and wherein said 30 
first successive spectral line comprises one of a 
sequence of spectral lines occurring during a sec- 
ond, contiguous time interval. 



6. The ultrasonic diagnostic apparatus of Claim 1, 35 
wherein said means for utilizing comprises 
means for identifying a velocity data value which 
is in the vicinity of the highest velocity intersec- 
tion of said spectral line data and said threshold 
level. 40 



7. The ultrasonic diagnostic apparatus of Claim 1, 
wherein said means for displaying comprise 
means for displaying a line connecting succes- 
sive identified peak velocity values. 45 

8. The ultrasonic diagnostic apparatus of Claim 2, 
wherein said means for displaying comprise 
means for displaying said spectral lines in a first 
color or tone, said identified peak velocity values so 
in a color or tone which is distinguished from said 

first color or tone, and said identified mean veloc- 
ity values in a color or tone which is distinguished 
from that of said spectral lines and said peak ve- 
locity values. 55 

9. The ultrasonic diagnostic apparatus of Claim 8. 
wherein said means for displaying comprise 
means for displaying said spectral lines in shades 
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